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Abstract. Slow light in photonic crystals and other periodic structures is 
associated with stationary points of the photonic dispersion relation, where the 
group velocity of light vanishes. We show that in certain cases, the vanishing 
group velocity is accompanied by the so-called frozen mode regime, when the 
incident light can be completely converted into the slow mode with huge diverging 
amplitude. The frozen mode regime is a qualitatively new wave phenomenon - 
it does not reduce to any known electromagnetic resonance. Formally, the frozen 
mode regime is not a resonance, in a sense that it is not particularly sensitive 
to the size and shape of the photonic crystal. The frozen mode regime is more 
robust and powerful, compared to any known slow-wave resonance. It has much 
higher tolerance to absorption and structural imperfections. 
PACS: 42.25.-p, 42.70.-a, 42.70.Qs 
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1. Introduction 



1.1. What is slow light? 

The velocity of light in transparent dispersive media is characterized by two different 
physical quantities - the phase velocity and the group velocity. The phase velocity of 
a travelling wave is defined as 

Vph = w/fc, (1) 

where lu is the frequency and k is the wave number. The group velocity of a travelling 
wave is defined as 

vg = duj/dk. (2) 

The group velocity of light is one of the most important electromagnetic characteristics 
of transparent media. It determines the speed of pulse propagation and is usually 
referred to simply as the speed of light propagation. With certain reservations, 
the group velocity also coincides with electromagnetic energy velocity in a wave 
packet [2 [5] . The definitions H]) and ^ apply not only to uniform transparent media 
and waveguides, but also to nonuniform spatially periodic structures, such as photonic 
crystals, arrays of coupled optical resonators, etc. In the case of a nonuniform periodic 
structure, k is the Bloch wave number defined in the Brillouin zone [T]. In a slow light 
case, the electromagnetic pulse propagates through the dispersive medium with the 
speed Vg <^ c, while the respective value of the phase velocity ^ is irrelevant. 

The term "slow light" has been widely used to describe a broad range of 
qualitatively different physical phenomena. There have been hundreds of publications 
in which the above term appears as a key word, while in many cases, the physical 
phenomena in question have little in common. Below, we just briefly outline the most 
important situations in which the concept of slow light can be invoked, and then we 
turn to the main subject of this paper - the frozen mode regime. 

Most of the "slow light" cases can be grouped into two major categories: 

- those where the low group velocity of light results from strong temporal 
(frequency) dispersion of the transparent medium; 

- those where the low speed of pulse propagation is a result of coherent interference 
in spatially periodic structures. 



1.2. Slow light in media with strong temporal dispersion 

In a uniform medium, both the phase and the group velocities of light can be expressed 
in terms of the refractive index n 

Vf^c/n, (3) 
+ ^d^) • 

The derivative dn/duj in ([4]) characterizes the frequency dispersion of the medium. 
At optical frequencies, the refractive index n of common transparent substances does 
not exceed several units, and the speed of light propagation is of the same order 
of magnitude as in vacuum. The situation can change dramatically in a strongly 
dispersive medium where the term ujdn/duj in ^ becomes dominant 

dn f dn\ ^ 

if Lo— ^ n, then w„ « c uj— ^ c. (5) 
dui \ duj J 
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Strong frequency dependence of the refractive index n implies that the hght group 
velocity can be substantially different from c. The relation ^ yields the following 
limitation on the slow pulse bandwidth 



where the assumption is made that the refractive index n within the transparency 
window is of the order of unity. The condition Q can also be viewed as a constraint 
on the minimal propagation speed of a pulse with a given bandwidth Aw. 

Usually, the large temporal dispersion ([5]) is a result of excitation of electronic 
or some other intrinsic resonances of the medium. This always involves some kind 
of a delicate resonant light-matter interaction with extremely small bandwidth. 
For instance, in the well known case of electromagnetically induced transparency 
m m III [3 O [TOj , the incident light interacts with atomic spin excitations forming 
the so-called dark-state polaritons. These polaritons propagate slowly through the 
medium in the form of a sharply compressed pulse, the energy of which is much 
smaller than that of the incident light pulse. In another example, the slow wave is 
associated with plasmons, which are electron density waves in gas of mobile electrons 
of the metal. Plasmons can also interact with light forming polaritons. There are 
many more examples of condensed matter excitations, strongly interacting with light 
and having low group velocity and relatively low relaxation rate. In most of such 
cases, the slow pulse propagating through the dispersive medium can be seen as an 
intrinsic coherent excitation (a polariton) triggered by the input light, rather than a 
light pulse per se. The group velocity of such an excitation may have little to do with 
the Maxwell equations and the speed of light c. 

In the rest of the paper we focus exclusively on the cases not involving any intrinsic 
excitations of the medium and, therefore, not related to temporal dispersion. Instead, 
we will focus on spatially periodic structures, in which low speed of electromagnetic 
pulse propagation results solely from spatial inhomogeniety. 

1.3. Slow light in spatially periodic structures 

Well-known examples of periodic dielectric structures include photonic crystals [HI [121 
[T1[T1], periodic arrays of coupled optical resonators [TCI [T7 1 [T5| . [1^1 . 
and photonic crystal coupled waveguide [2l[25l[26l[2Zl[28l[29l[30l[3T]. Low speed 
of pulse propagation in all these cases is the result of coherent interference of light 
scattered at the interfaces of adjacent structural components. The effects of spatial 
dispersion associated with structural periodicity are particularly strong when the 
structural period L and the light wavelength A are comparable in value 



Strong spatial dispersion can result in low or even zero group velocity of the 
respective Bloch wave. In practice, at optical frequencies, the speed of pulse 
propagation in periodic dielectric arrays can be reduced by not more than three 
orders of magnitude. This is not a fundamental restriction, but rather a technological 
limitation related to the difhculty of building flawless periodic arrays at nanoscales. 
On the positive side, the dielectric components of the periodic array are not required 
to display strong temporal dispersion. As a consequence, the absorption of light is not 
an inherent problem in this case. 




(6) 



A. 



(7) 
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There is a natural bandwidth limitation on the slow pulse in periodic dielectric 
arrays, which is similar to the case of slow light in time-dispersive media. Let Aw be 
the frequency bandwidth of a pulse and AA; - the respective range of the Bloch wave 
number. The average group velocity (vg) of the pulse can be defined as 

Let us assume that the pulse propagating in the periodic medium is composed of the 
Bloch eigenmodes belonging to the same spectral branch of the dispersion relation 
Lo (k). This assumption implies that Afc cannot exceed the size 2t:/L of the Brillouin 
zone 

Afc < 27r/i, 

where L is the unit cell length of the periodic array. In addition, we assume that the 
refractive index of the constitutive components of the periodic array is of the order of 
unity and, therefore, 

L - Ao = 27rc/w, (9) 

where Aq is the light wavelength in vacuum. The relations ([5]) and ^ yield the 
following limitation on the minimal propagation speed of a pulse with a given 
bandwidth Ao; 

{v^)>IlAu;^c—. (10) 

ZTT LU 

The restriction (jTU]) is similar to that defined by the inequality ^ and related to the 
case of slow light in a uniform medium with temporal (frequency) dispersion. 



1.4- Examples of periodic arrays supporting slow light 

During the last two decades, a tremendous progress has been made in theory and 
apphcations of periodic arrays of coupled optical resonators. Generally, if the coupling 
between adjacent resonators in a periodic array is weak, the group velocity of Bloch 
excitations propagating through the array is very low in each and every transmission 
band. This is true regardless of the nature of individual resonators. The above simple 
idea forms the basis for one of the most popular approaches to slowing down the light. 
One inevitable consequence of the weak coupling between the neighboring optical 
resonators is that all the individual transmission bands are very narrow. Another 
problem is that there can be severe restrictions on transmitted power. An extensive 
discussion on the subject and numerous examples and references can be found in 
[igilllllllllHllIglEQlETlEllEa]. a qualitatively similar situation occurs in line- 
defect waveguides in photonic crystals, where a periodic array of structural defects 
plays the role of weakly coupled optical resonators [M l [?7 l 1 ^ 1 ^ I^H l [5T] . 

Slow light phenomena in periodic arrays of weakly coupled resonators have been 
the subject of a great number of recent publications, including some excellent review 
articles cited above. A common characteristic of all different realizations of this 
approach is a relatively low density of modes. One consequence of this is a significant 
nonlinearity. In certain cases, the nonlinearity is so extreme that it occurs on a single 
photon level. Strong nonlinearity can be useful in controlling the flow of light. On the 
other hand, it can severely limit the transmission capabilities of the optical waveguide. 
In this respect, photonic crystals have advantage. 
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Photonic crystals are spatially periodic structures composed of two or more 
different transparent dielectric materials [TTl [T^l [131 |T3] . Unlike the case of optical 
waveguides and linear arrays of coupled resonators, in photonic crystals we have bulk 
electromagnetic waves capable of propagating in any direction through the periodic 
heterogeneous structure. This results in much greater density of modes, compared to 
arrays of weakly coupled resonators [32l [33l [34j . 

Strong spatial dispersion of a typical photonic crystal and any other periodic 
dielectric structure is reflected in a complicated k — uj diagram featuring transmission 
bands and gaps. Normally, each spectral branch ui {k) of the Bloch dispersion 
relation develops stationary points uis — (ks) where the group velocity ([2]) of the 
corresponding propagating mode vanishes 

duj/dk ~0, atuj — LUs—uj (kg) . (11) 

Examples of different stationary points are shown in Fig. [TJ Each of the frequencies tOg, 
ivo and lu^ is associated with zero group velocity of the respective traveling wave, but 
there are some fundamental differences between those three cases. These differences 
are particularly pronounced when it comes to the efficiency of conversion of the 
incident light into the slow mode inside the periodic structure. In most cases, incident 
radiation with the frequency of one of the slow modes is simply reflected back to 
space, without exciting the slow mode inside the periodic structure. How to overcome 
this fundamental problem and, thereby, how to transform a significant fraction of the 
incident light energy into a slow mode with drastically enhanced amplitude, is the 
primary subject of this paper. 

2. Frozen mode regime in photonic crystals. 

The effects of strong spatial dispersion in periodic dielectric structures culminate at 
stationary points (fTT|) of the Bloch dispersion relation, where the group velocity ([2]) of 
a traveling Bloch wave vanishes. One reason for this is that vanishing group velocity 
always implies a dramatic increase in density of modes at the respective frequency. 
In addition, vanishing group velocity also implies certain qualitative changes in the 
eigenmode structure, which can strongly affect the propagation and scattering of 
light. A particular example of the kind is the frozen mode regime associated with 
a dramatic enhancement of the amplitude of the wave transmitted to the periodic 
medium [35, 36, 37, 38 . There are at least two qualitatively different modifications of 
the frozen mode regime, each related to a specific singularity of the electromagnetic 
dispersion relation. Both effects can be explained using a simple example of a plane 
electromagnetic wave normally incident on a lossless semi-infinite periodic structure. 

The frozen mode regime of the first kind is associated with a stationary infiection 
point (SIP) of the k — u diagram shown in Fig. [ijb). In the vicinity of stationary 
inflection point, the relation between the frequency u and the Bloch wave number k 
can be approximated as 

Lu — luq (X {k ~ ko)'^ . (12) 

A monochromatic plane wave of frequency close to loq incident on semi-infinite 
photonic crystal is converted into the frozen mode with infinitesimal group velocity 
and dramatically enhanced amplitude, as illustrated in Fig. [21 The saturation value of 
the frozen mode amplitude diverges as the frequency approaches the SIP value uq, as 
shown in Fig. [Ifc). Remarkably, the photonic crystal reflectivity at w = can be very 
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low, implying that the incident radiation is almost totally converted into the frozen 
mode with zero group velocity, drastically enhanced amplitude, and finite energy flux 
close to that of the incident wave. This remarkable phenomenon is uniquely associated 
with a stationary inflection point (SIP) of the Bloch dispersion relation. Not every 
periodic structure can display such a spectral singularity, but if it does, one can expect 
the behavior described in Fig. [21 

A qualitatively different kind of frozen mode regime can occur in the vicinity of a 
degenerate photonic band edge (DBE) shown in Fig. [He). At frequencies just below 
LOd, the respective dispersion relation can be approximated as 

Wrf — W CX (fc — fej;)"* , at W ^ LLJd- (13) 

Contrary to the case of stationary inflection point (|12p . in the vicinity of a degenerate 
band edge the photonic crystal becomes totally reflective. But at the same time, 
the steady-state field inside the periodic medium (at z > 0) develops a very large 
amplitude, diverging as the frequency approaches its critical value LOd- Such a 
phenomenon is illustrated in Fig. [3) The frozen mode profile below and above the 
degenerate band edge frequency ujd is different. It has a large saturation value at 
frequencies located inside the transmission band (at lu < ujd)i as seen in Fig. [3Ka) and 
(b). This saturation value diverges as lo —>■ ujd ~ 0. By contrast, at frequencies lying 
inside the band gap (at u > w^), the field amplitude initially increases dramatically 
with the distance z from the surface, but then vanishes as the distance z further 
increases, as seen in Fig. - f). 

Hereinafter, the above two phenomena will be referred to as the SIP-related frozen 
mode regime and the DBE-related frozen mode regime, respectively. Figs. [2] and [3] 
describe the frozen mode regime in hypothetical lossless semi-infinite photonic crystal. 
What happens to the frozen mode regime in a plane-parallel photonic slab with a finite 
thickness D? 

First, let us assume that the incident light frequency is equal to that of the 
respective stationary point {too or ujd) of the Bloch dispersion relation in Fig. [T] 
In either case, in the leftmost portion of the photonic slab, the frozen mode profile 
remains the same as in the semi- infinite case shown in Figs. HJc) and[3Jc). In the 
opposite, rightmost part of the photonic slab, the frozen mode amplitude now vanishes, 
as shown in Fig. |4Ka-c). In either case, the maximum field intensity Wmax inside the 
photonic slab is 

VFmax CC WlN^, (14) 

where Wj is the intensity of the incident light and N — D/L is the number of periods 
(primitive translations) in the z direction of the finite periodic structure. Additional 
factors limiting the frozen mode amplitude include structural imperfections of the 
periodic array, absorption, nonlinearity, deviation of the incident radiation from plane 
monochromatic wave, etc. 

Now, what happens in a finite photonic slab if the incident light frequency is 
slightly different from that of the respective stationary point loq or cj^? 

In the case of a SIP-related frozen mode regime, the field amplitude will decrease 
as the frequency lo deviates from loq in either direction ~ similar to the case of semi- 
infinite structure shown in Fig. [2l In other words, in a lossless finite photonic slab, 
the SIP-related frozen mode regime appears as a resonance centered at wq- According 
to (fH)) . the light intensity at the resonance is proportional to N'^, while the respective 
Q-factor is proportional to N^. 
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In the case of a DBE-related frozen mode regime, the situation is quite different. 
Namely, in a lossless, plane-parallel photonic sample with a finite thickness, the DBE- 
related frozen mode regime is overwhelmed with much more powerful, giant slow wave 
resonance [321 HO]- Let us explain this phenomenon in more detail. At DBE frequency 
u)d, the frozen mode profile in a finite photonic slab is shown in Fig. IHa-c). According 
to ([T4|) . the maximum light intensity inside the photonic slab is proportional to N'^. 
If the incident light frequency deviates from cod toward the photonic band gap in 
Fig. IHc), the light intensity inside the photonic slab will decrease, similar to the 
semi- infinite case shown in Figs. [3{d-f). But if the incident light frequency deviates 
from LUd in the opposite direction (toward the transmission band in Fig. [He)), the 
light intensity inside the photonic slab will increase dramatically, peaking at a certain 
resonance frequency ujr, as shown in Fig. [5l The exact location of the resonance 
frequency tOr is dependent on the thickness of the finite periodic structure |39| |40| . But 
in any event, ujr lies inside the transmission band and approaches the DBE frequency 
LOd as N increases. The light intensity at the frequency uJr of the giant slow wave 
resonance is 

VKmax K WlN\ (15) 

which is by factor iV^ larger than that of the common slow wave resonance associated 
with the regular photonic band edge in Fig. [BJa). The Q-factor of the giant slow wave 
resonance is proportional to N^. 

In summary, in a photonic slab of finite thickness, the SIP-related frozen mode 
regime persists, but its amplitude is limited by the slab thickness in accordance with 
the relation (fT4|) . By contrast, the DBE-related frozen mode regime in a finite plane- 
parallel photonic slab is overrun by the giant slow wave resonance, which can be 
viewed as a combination of the frozen mode regime and the Fabry-Perot resonance. 
It is much more powerful compared to the common RBE-related slow wave resonance 
in a photonic slab of similar size [39l |40] . 

Not every periodic structure can support the frozen mode regime at normal 
incidence. Generally, the physical conditions for the frozen mode regime are the 
same as the conditions for the existence of the respective stationary points (fT^ or 
(|13p of the Bloch dispersion relation. For instance, in the case of periodic layered 
structure, a unit cell must contain at least three layers, of which two must display 
a misaligned in-plane anisotropy [45l [35l [38] . At optical frequencies, fabrication of 
such a layered structure can be challenging. Fortunately, in photonic crystals with 3D 
periodicity, the conditions for the existence of stationary inflection point (SIP) and 
the related frozen mode regime are much less restrictive. A good example of the kind 
is provided by the inverse opal structure, where such a spectral singularity was found 
and investigated, both theoretically and experimentally in [46| . 

The frozen mode regime is a qualitatively new wave phenomenon - it does not 
reduce to any known electromagnetic resonance. Formally, the SIP-related frozen 
mode regime is not a resonance, in a sense that it is not particularly sensitive to 
the size and shape of the photonic crystal. Besides, the frequency dependence of the 
frozen mode amplitude is very different from that of a cavity resonance, or a common 
slow wave resonance. The frozen mode regime is much more robust, compared to any 
known slow-wave resonance occurring in periodic and nonperiodic photonic structures. 
It has much higher tolerance to absorption and structural imperfections than common 
Fabry-Perot or transmission band edge resonances, where the entire photonic crystal 
works as a resonator. 
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The DBE-related frozen mode regime is quite different. In a finite plane-parallel 
slab it is overrun with the much more powerful giant slow- wave resonance [311 HHI- The 
Q-factor associated with such a resonance can be by two orders of magnitude higher, 
compared to that of the regular slow-wave resonance in the same or similar periodic 
structure. More importantly, it provides the possibility of a dramatic reduction in size 
- up to an order of magnitude - of some basic photonic devices without compromising 
their performance. 



3. Physical nature of the frozen mode regime 



The essence of the frozen mode regime can be understood from the simple example 
of a plane monochromatic wave normally incident on a semi-infinite periodic layered 
structure. An important requirement, though, is that some of the layers display a 
misaligned in-plane anisotropy [151 [351 13H] ■ Below we present a comparative analysis 
of two different kinds of frozen mode regime. Although throughout this section we 
only consider the case of normal incidence, most of the results and expressions remain 
virtually unchanged in a more general case of the frozen mode regime at oblique 
propagation [36l [3ll [38]. One difference, though, is that at oblique incidence, the 
frozen mode regime can occur in much simpler structures, which can have a big 
practical advantage. 

To start with, let us introduce some basic notations and definitions. Let ^f/, ^fjj, 
and \E't be the incident, reflected and transmitted waves, respectively. We assume 
that all three monochromatic waves propagate along the z axis normal to the layers. 
Electromagnetic field both inside (at z > 0) and outside (at z < 0) the semi-infinite 
periodic structure is independent of the x and y coordinates. The transverse field 
components can be represented as a column-vector 



vl'(z) 



Exiz) 

Ey [Z) 

Hy{z) 



(16) 



where E (z) and H (z) are time-harmonic electric and magnetic fields. All four 
transverse field components in are continuous functions of z, which leads to the 
following standard boundary condition at the photonic crystal interface at z = 

*T (0) - */ (0) + (0) . (17) 

Assume also that anisotropic layers of the periodic array have an in-plane birefringence 

^xx ^x 



-xy 





-xy 

-yy 







^ zz 



in which case the fields E (z) and H (z) are normal to the direction of light propagation 

E{z) ±z,H{z) ±z. (18) 



Note that the polarizations of the incident, reflected, and transmitted waves can be 
all different, because some of the layers of the periodic array are birefringent. The 
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presence of birefringent layers is essential for the possibility of frozen mode regime. 
Moreover, in the case of normal incidence, each unit cell of the periodic stack must 
include at least two birefringent layers with misaligned anisotropy axes [SS] [321 [3H] ■ 

In periodic layered media, the electromagnetic eigenmodes (z) are usually 
chosen in the Bloch form 

*fc(z + L) = e*'^'^*fe(z), (19) 

where the Bloch wave number k is defined up to a multiple of 2-k/L. The 
correspondence between lo and k is referred to as the Bloch dispersion relation. Real 
wave numbers k correspond to propagating (traveling) Bloch modes. Propagating 
modes belong to different spectral branches w [k) separated by frequency gaps. The 
speed of a traveling wave in a periodic medium is determined by the group velocity 
([2]). Complex wave numbers k = k' + ik" correspond to evanescent Bloch modes. 
Evanescent modes decay exponentially with the distance z from the boundary of 
semi-infinite periodic structure. Therefore, under normal circumstances, evanescent 
contribution to the transmitted wave ^'t {z) can be significant only in close proximity 
of the photonic crystal interface at z = 0. The situation can change dramatically 
when the frequency a; approaches one of the stationary point values lOs- At first 
sight, stationary points (fTT|) relate only to propagating Bloch modes. But in fact, 
in the vicinity of every stationary point frequency Wg, the imaginary part k" of the 
Bloch wave number of at least one of the evanescent modes also vanishes. As a 
consequence, the respective evanescent mode decays very slowly, and its role may 
extend far beyond the photonic crystal boundary. In addition, in the special cases of 
interest, the electromagnetic field distribution ^E* (z) in the coexisting evanescent and 
propagating eigenmodes becomes very similar, as uj approaches Us- This can result in 
an abnormal interference pattern constituting the frozen mode regime. What exactly 
happens in the vicinity of a particular stationary point (jlll) essentially depends on its 
character and appears to be very different in each of the three cases presented in Fig. 

m 

3.1. Bloch composition of frozen mode 

If all the layers of the periodic layered structure were made of isotropic materials, such 
as glass or air, each of the Bloch eigenmodes would be doubly degenerate with respect 
to light polarization. This applies both to propagating and evanescent Bloch modes. 
In our case, though, each unit cell of the periodic array contains a pair of birefringent 
layers with misaligned anisotropy axes. As a consequence, the Bloch eigenmodes 
with different (elliptic) polarization will not be degenerate. At any given frequency, 
the total number of Bloch eigenmodes with k \\ z is four. The only exception are 
the stationary point frequencies (|1H) . where some of the four eigenmodes cannot be 
represented in the Bloch form (fT9|) . Those cases will be discussed later. 

With the exception of stationary points (jlll) , any plane monochromatic wave ()16|) 
can be represented as a superposition of four Bloch eigenmodes, propagating and/or 
evanescent, with different polarizations and wave numbers. But in the setting where 
the semi- infinite periodic array occupies the half-space z > 0, the transmitted wave 
is a superposition of only two of the four Bloch eigenmodes. Indeed, neither the 
propagating modes with negative group velocity, nor evanescent modes exponentially 
growing with the distance z from the surface, contribute to ^'t (z) in this case. 
Generally, one can distinguish three different possibilities. 
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(i) Both Bloch components of the transmitted wave ^'^ are propagating modes 

*T (^) = *prl (-2) + *P'-2 (2) , Z>0. (20) 

^pri (z) and ^pr2 (z) are two propagating Bloch modes with different real wave 
numbers ki and fc2 and different group velocities Vgi > and Vg2 > 0. This 
constitutes the phenomenon of double refraction, provided that Vgi and Vg2 are 
different. The remaining two Bloch modes of the same frequency have negative 
group velocities and cannot contribute to the transmitted wave \1/t- 

(ii) Both Bloch components of '^t are evanescent 

^T{z) = -^evl{z) + ^ev2{z), Z>0. (21) 

The respective two values of k are complex with positive imaginary parts k'l > 0, 
^2 > 0. This is the case when the frequency ui falls into photonic band gap at 
uj > ujg in Fig. [DJa) or at w > U!d in Fig. [TKc). The fact that k'{ > 0, > 
implies that the wave amplitude decays with the distance z from the surface. In 
the case ([2T|). the incident wave is totally reflected back to space by the semi- 
infinite periodic structure. 

(iii) One of the Bloch components of the transmitted wave ^'t is a propagating mode 
with Vg > 0, while the other is an evanescent mode with k" > 

(z) = "^fpr (z) + ^ev (z), Z> 0. (22) 

For example, this is the case a.t oj ^ ujq in Fig. [Hb), as well as at w < ujg in Fig. 
[ija) and at uj < Wd in Fig. [TJc). As the distance z from the surface increases, 
the evanescent contribution \l/e^ in ([^ decays as exp {—zk"), and the resulting 
transmitted wave '^t (z) turns into a single propagating Bloch mode "^pr- 

Propagating modes with Vg > and evanescent modes with k" > are referred 
to as forward modes. Only forward modes contribute to the transmitted wave ^'t (z) 
in the case of a periodic semi-infinite stack. The propagating modes with Vg < and 
evanescent modes with k" < are referred to as backward modes. The backward 
Bloch modes never contribute to the transmitted wave (z) inside the periodic 
semi-infinite stack. This statement is based on the following two assumptions: 

- The transmitted wave and the reflected wave ^'^ are originated from the 
plane wave "ifj incident on the semi-infinite photonic slab from the left. 

- The periodic structure occupies the entire half-space and is perfectly periodic at 
2 > 0. 

If either of the above conditions is violated, the electromagnetic field inside the 
periodic stack can be a superposition of four Bloch eigenmodes with either sign of 
the group velocity Vg of propagating contributions, or either sign of k" of evanescent 
contributions. This would be the case if the periodic array had some kind of structural 
defects or a finite thickness. At the end of this section we briefiy discuss how it would 
affect the frozen mode regime. 

Note also that the assumption that the transmitted wave ^'t (z) is a superposition 
of propagating and/or evanescent Bloch eigenmodes may not apply if the frequency 
Lv exactly coincides with one of the stationary point frequencies (jlip . For example, 
at frequency ujq of stationary inflection point (fT2)) . there are no evanescent solutions 
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to the Maxwell equations, and the transmitted wave (z) is a (non-Bloch) Floquet 
eigenmode linearly growing with z [321 [M| ■ Similar situation occurs at frequency cud 
of degenerate band edge (fT5)) . The term "non-Bloch" means that the respective field 
distribution does not comply with the relation (|19p . At the same time, at any general 
frequency, including the vicinity of any stationary point (jll[) . the transmitted wave 
(z) is a superposition of two forward Bloch eigenmodes, each of which is either 
propagating, or evanescent. 

In all three cases (f20l - [22]) , the contribution of a particular Bloch eigenmode to 
the transmitted wave ^'t depends on the polarization of the incident wave. One 
can always choose some special (elliptic) incident wave polarization, such that only one 
of the two forward Bloch modes is excited and the transmitted wave is a single 
Bloch eigenmode. In the next subsection we will see that there is no frozen mode 
regime in the case of a single mode excitation. This fact relates to the very nature of 
the frozen mode regime in a periodic layered structure [38] . 

Knowing the Bloch composition of the transmitted wave we can give a semi- 
qualitative description of what happens when the frequency to of the incident wave 
approaches one of the stationary points (fTT|) in Fig. [TJ The rigorous mathematical 
analysis of the respective scattering problem can be found in [35l [38] . 

3.1.1. Regular photonic band edge We start with the simplest case of a regular 
photonic band edge. There are two different possibilities in this case, but none of 
them is associated with the frozen mode regime. The first one relates to the trivial 
case where none of the layers of the periodic structure displays an in-plane anisotropy 
or gyrotropy. In this case, all eigenmodes are doubly degenerate with respect to 
polarization. A detailed description of this case can be found in the extensive literature 
on optics of stratified media [121 [43] . Slightly different scenario occurs if some of the 
layers are anisotropic or gyrotropic and, as a result, the polarization degeneracy is 
lifted. Just below the band edge frequency LOg in Fig. [TJa), the transmitted field 
[z) is a superposition ([22]) of one propagating and one evanescent Bloch modes. 
Due to the boundary condition (jl7[) , the amplitude of the transmitted wave at z = 
is comparable to that of the incident wave. In the case of a generic polarization of the 
incident light, the amplitudes of the propagating and evanescent Bloch components 
at z = are also comparable to each other and to the amplitude of the incident light 



As the distance z from the surface increases, the evanescent component ev {z) 
decays rapidly, while the amplitude of the propagating component remains constant. 
Eventually, at a certain distance from the slab surface, the transmitted wave ^'t [z) 
becomes very close to the propagating mode 



The evanescent component of the transmitted wave does not display any 
singularity at the band edge frequency ujg. The propagating mode ^'p^ does develop a 
singularity associated with vanishing group velocity at a; o^g — 0, but its amplitude 
remains finite and comparable to that of the incident wave. At w > Wg, this 
propagating mode turns into another evanescent mode in (|2ip . The bottom line is that 
none of the Bloch components of the transmitted wave develops a large amplitude in 
the vicinity of a regular photonic band edge. There is no frozen mode regime in this 
trivial case. 



\^pr (0)1 - |*e. (0)1 - |*/| , at < Ug. 



(23) 



^-T (z) « ^'pr (z) , at Z > i, UJ < U)g. 



(24) 
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3.1.2. Stationary inflection point A completely different situation develops in the 
vicinity of a stationary inflection point (HH). At a; « cjqi the transmitted wave 

is a superposition (j22p of one propagating and one evanescent Bloch component. 
In contrast to the case of a regular photonic band edge, in the vicinity of loq both 
Bloch contributions to develop strong singularity. Specifically, as the frequency 
uj approaches , both contributions grow dramatically, while remaining nearly equal 
and opposite in sign at the slab boundary [35j 

^fpr (0) « -^fet, (0) oc |w - woT^^^ , a.s uj^ujo. (25) 
Due to the destructive interference the resulting field 

*T (0) = V (0) + ^e. (0) 

at the surface at z = is small enough to satisfy the boundary condition (|17p . As the 
distance z from the slab boundary increases, the destructive interference becomes less 
effective - in part because the evanescent contribution decays exponentially 

■^ev [z) « -^ev (0) exp {-zk") , (26) 

while the amplitude of the propagating contribution remains constant and very 
large. Eventually, the transmitted wave V^t (z) reaches its large saturation value 
corresponding to its propagating component ^'pr, as seen in Fig. [6l^a). 

Note that the imaginary part k" of the evanescent mode wave number in (|26|1 
also vanishes in the vicinity of stationary inflection point 

k" oc \lu — luqI^^^ , as w ^ wq; (27) 

reducing the rate of decay of the evanescent contribution (|26|) . As a consequence, the 
resulting amplitude ^'t (z) of the transmitted wave reaches its large saturation value 
\E'pr in (|25p only at a certain distance Z from the surface 

Zocl/fc" oc |w-wo^^/^ (28) 

This characteristic distance diverges as the frequency approaches its critical value loq . 

If the frequency of the incident wave is exactly equal to the frozen mode frequency 
Wo, the transmitted wave ^'t {z) does not reduce to the sum (|22|) of propagating and 
evanescent contributions, because at lo — loq, there is no evanescent solutions to the 
Maxwell equations. Instead, \1/t (z) corresponds to a non-Bloch Floquet eigenmode 
diverging linearly with z |35) . 

*T (z) - *T (0) oc z^o, atw=a;o- (29) 
Such a solution is shown in Fig. [SJc). 

3.1.3. Degenerate band edge While the situation with a regular photonic band edge 
appears trivial, the case of a degenerate band edge ifTS]) proves to be quite different. 
Just below the degenerate band edge frequency lu^ (inside the transmission band), 
the transmitted field is a superposition ([22l) of one propagating and one evanescent 
components. Above Ud (inside the band gap), the transmitted wave is a combination 
(|2ip of two evanescent components. In this respect, a regular and a degenerate band 
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edges are similar to each other. A crucial difference, though, is that in the vicinity 
of a degenerate band edge, both Bloch contributions to the transmitted wave diverge 
as Lu approaches both above and below the band edge frequency. This constitutes 
the frozen mode regime. 

Let us start with the transmission band. As the frequency uj approaches Ud — 0, 
both Bloch contributions in (|22|) grow sharply, while remaining nearly equal and 
opposite in sign at the surface at z = [38] 

*p^(0) W -*et,(0) OC a.S LU (30) 

The destructive interference ((30|) ensures that the boundary condition (fT7| can be 
satisfied, while both Bloch contributions to (z) diverge. As the distance z from 
the slab boundary increases, the evanescent component \E'e« (z) dies out 

*e«(z)«*e.(0)exp(-zfc") (31) 

while the propagating component '^pr (z) remains constant and very large. Eventually, 
as the distance z further increases, the transmitted wave {z) reaches its large 
saturation value corresponding to its propagating component (z), as illustrated in 
Fig. [T] Note that the imaginary part k" of the evanescent mode wave number also 
vanishes in the vicinity of degenerate band edge 

k" oc — ujdl'^^'^ , as u —I- Ud, (32) 

reducing the rate of decay of the evanescent contribution (|3ip . As a consequence, the 
resulting amplitude ^'t (z) of the transmitted wave reaches its large saturation value 
^pr only at a certain distance Z from the surface 

Z ocl/k" oc\uj-uJd\^^^^ ■ (33) 

This characteristic distance increases as the frequency approaches its critical value LOd, 
as illustrated in Fig. [3Ka) and (b). 

If the frequency w of the incident wave is exactly equal to uid, the transmitted wave 
\1/t {z) does not reduce to the sum of two Bloch contributions. Instead, it corresponds 
to a non-Bloch Floquet cigcnmode linearly diverging with z 

{z) - *T (0) OC z^d, at a; = ujd- (34) 

This situation is shown in Fig. [3Kc). 

The above behavior appears to be very similar to that of the frozen mode regime 
at a stationary inflection point, shown in Figs. [2]and[6l Yet, there is a crucial difference 
between the frozen mode regime at a stationary inflection point and at a degenerate 
band edge. In the immediate proximity of a degenerate band edge, the Pointing vector 
5*^ of the transmitted wave is infinitesimal, in spite of the diverging wave amplitude. 
In other words, although the energy density Wt oc \'^t\^ of the frozen mode diverges 
as a; — > — 0, it does not grow fast enough to offset the vanishing group velocity. As 
a consequence, the photonic crystal becomes totally reflective at uj — ujd- Of course, 
the total reflectivity persists at a; > ujd, where there is no propagating modes at all. 
By contrast, in the case (p9)) of a stationary inflection point, the respective Pointing 
vector St is finite and can be even close to that of the incident wave, implying low 
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reflectivity and nearly total conversion of the incident wave energy into the frozen 
mode. 

The character of frozen mode regime is different when we approach the degenerate 
band edge frequency from the band gap. In such a case, the transmitted field (z) 
is a superposition ([2T|l of two evanescent components. As the frequency uj approaches 
LVd, both evanescent contributions grow sharply, while remaining nearly equal and 
opposite in sign at the photonic crystal boundary 

^'eul (0) « -*e«2 (0) CX |cJd -t^r^^\ aS W ^ Wrf + 0. (35) 

Again, the destructive interference ([55]) ensures that the boundary condition ([TT]) can 
be satisfied, while both evanescent contributions to 'I't (z) diverge in accordance with 
(|35p . As the distance z from the slab boundary increases, the destructive interference 
of these two evanescent components is lifted and the resulting field amplitude increases 
sharply, as seen in Fig. [H^a). But eventually, as the distance z further increases, 
the transmitted wave (z) completely decays, because both Bloch contributions to 
(z) are evanescent. The latter constitutes the major difference between the frozen 
mode regime above and below the DBE frequency Ud- The rate of the amplitude 
decay, as well as the position of the maximum of the transmitted wave amplitude in 
Figs. [Sljd - f) and[8l^a) , are characterized by the distance Z in ((33)l . 

3.1.4- Physical reason for the growing wave amplitude If the frequency lu is close, 
but not equal, to that of a stationary point (jlip of the dispersion relation, the wave 
V&T (z) transmitted to the semi-infinite periodic layered structure is a superposition of 
two forward Bloch modes 5*1 (z) and (z) 

(z) = *i (z) + *2 (z) , z> 0. (36) 

The two Bloch modes in ((36|) can be one propagating and one evanescent, as in ((22)) . 
or they can be both evanescent, as in (|2ip . In the vicinity of frozen mode regime, as 
the frequency approaches its critical value of ujo or LOd, the four-dimensional vectors 
(|16p corresponding to each of the two Bloch eigenmodes on the right-hand side of ([36|) 
become nearly parallel to each other 

*i (z) « a*2 (2) , asw-^ws, (37) 

where a is a scalar, and cOg is the frozen mode frequency ujq or ujd. The asymptotic 
relation ([57]) reflects a basic property of the transfer matrix of the periodic layered 
structure at the frequency of either stationary point ujq or uJd (see [36-38]). 

Let us show how the property (|37p of the Bloch eigenmodes in the vicinity of ujq 
or UJd can lead to the frozen mode regime. Indeed, at the photonic crystal boundary, 
the sum p6|) of two nearly parallel column vectors ^'i (z) and 4*2 (z) must match the 
boundary conditions ([T7]) with the incident and reflected waves ^/ (0) and 'f/j (0). If 
the incident wave polarization is general, we have no reason to expect the column- 
vector 

*(0) = */(0) + *h(0) (38) 

at the photonic crystal interface to be nearly parallel to 'i'l (0) and ^'2 (0). But on 
the other hand, the boundary condition pT|) says that 



^^l (0) = 5-1 (0) + ^-2 (0) 



(39) 
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Obviously, the only situation where 5* (0) in (|38p and ([55)1 is a general four-dimensional 
vector, while ^'i (0) and (0) on the right-hand side of (|39p are nearly parallel to 
each other is when 

*i(0)«-*2(0), 1*1 (0)1 « 1*2(0)1 > I* (0)1 . (40) 

The relation (j40[) explains how the huge diverging amplitude of the frozen mode inside 
the photonic crystal can be reconciled with the boundary conditions (|17p . 

There is one exception, though. As we already stated in (I37p . in the vicinity 
of the frozen mode frequency, the two Bloch components *i (z) and *2 (2) of the 
transmitted wave *t (2) are nearly parallel to each other. For this reason, if the 
polarization of the incident wave */ is such that * (0) in p9|) is nearly parallel to 
one of the Bloch eigenmodes *i (0) or *2 (0), it is also nearly parallel to the other 
one. So, all three column vectors *i (0), *2 (0), and * (0) are now parallel to each 
other. In this, and only this case, the amplitude of the transmitted wave *t (z) will 
be comparable to that of the incident wave. There is no frozen mode regime for the 
respective vanishingly small range of the incident wave polarization. A particular case 
of the above situation is the regime of a single mode excitation, where only one of the 
two Bloch components *i (z) or *2 (z) in ([55]) contributes to the transmitted wave 
(z) [3E]- Indeed, if the transmitted wave is a single Bloch mode (either *i (z), or 
*2 (z)), the condition cannot be satisfied. 

Finally, let us reiterate that in the limiting cases of uj ~ loq or lo — ujd-, the 
transmitted wave *t (2) corresponds to the non-Bloch Floquet eigenmode ((29)l or 
P4[). respectively. Either of them linearly diverges with z. The only exception is 
when the incident wave has the unique (elliptic) polarization, for which the transmitted 
wave *T (z) is a propagating Bloch eigenmode with zero group velocity and a limited 
amplitude, comparable to that of the incident wave. Incident wave with any other 
polarization will generate the frozen mode inside the periodic medium. 

4. Summary 

Although the existence of slow electromagnetic modes in photonic crystals is quite 
obvious, the next question is whether and how such modes can be excited by incident 
light. We have shown that the frozen mode regime provides unique advantages in this 
respect. Generally, the possibility of the frozen mode regime is determined by the 
character of the Bloch dispersion relation, rather than the specific physical realization 
of the periodic structure supporting such a dispersion relation. As soon as the Bloch 
dispersion relation displays the proper singularity such as a stationary inflection point 
(|12l) or a degenerate band edge (|13p . we have every reason to expect the occurrence 
of the frozen mode regime at the respective frequency. In other words, the possibility 
of the frozen mode regime is determined by the spectral properties of the periodic 
structure, rather than by specific physical nature of the linear excitations. 

If a periodic structure is relatively simple - for instance, a stratified medium with 
one dimensional periodicity - its frequency spectrum may prove to be too simple to 
support the proper spectral singularity and the frozen mode regime. All our numerical 
examples relate to ID photonic crystals, where the existence of the frozen mode regime 
requires the presence of misaligned birefringent layers (see the details in 3F). The 
more complex the periodic structure is, the more likely it is capable of supporting such 
a phenomenon. In photonic crystals with 3D periodicity, there are no fundamental 
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restrictions on the existence of the frozen mode regime |46) . The same is true for the 
cases of modulated waveguides and periodic arrays of coupled resonators. 

Another important question is how robust the frozen mode regime is. For 
instance, what happens if we introduce absorption or structural imperfections. Of 
course, all these factors suppress the frozen mode amplitude, but not as much as in 
the case of common Fabry-Perot or transmission band edge resonances, where the 
entire photonic structure works as a resonator. In addition, since in the vicinity of an 
inflection point the dispersion term v'g (k) — lo" (fc) vanishes, one can achieve wide- 
bandwidth and dispersion-free propagation of light [27] . 

Finally, it is possible to combine the frozen mode regime and a common slow 
wave resonance, in which case we have the phenomenon called the giant transmission 
resonance [39l |40]. The Q-factor associated with such a resonance can be by two 
orders of magnitude higher, compared to that of the regular Fabry-Perot resonance in 
the same or similar periodic structure. 
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Figure 1. (Color online). Schematic examples of dispersion relations displaying 
different stationary points: (a) a regular band edge (RBE), (b) a stationary 
inflection point (SIP), (c) a degenerate band edge (DBE). 
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Figure 2. (Color online). Smoothed profile of the frozen mode at six different 
frequencies in the vicinity of stationary inflection point: (a) uj = ljq — W~'^c/L, 
(b) Lu = uJo~ IQ-^c/L, (c) uj = Luo, (d) lu = uJo + lQ~^c/L, (e) uj = cJo + 10"*c/L, 
(f) UJ = ujQ + 10~'^ c/L. In all cases, the incident wave has the same polarization 
and unity amplitude. The distance z from the surface of semi-infinite photonic 
crystal is expressed in units of L. 
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Figure 3. (Color online). Smoothed profile of the frozen mode at six different 
frequencies in the vicinity of degenerate band edge: (a) uj = u>d — 10~*c/L, (b) 
w = - lO-^c/L, (c) Ld = uJd, (d) uj = LUa + 10-^c/L, (e) w = uja + IQ-^c/L, 
(f) u! = + 10~'^c/L. In the transmission band (at uj < iOd)y the asymptotic 
field value diverges as a; — > ujj. By contrast, in the band gap (at uj > uid), the 
asymptotic field value is zero. The amplitude of the incident wave at z < is 
unity. The distance z from the surface is expressed in units of L. 
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Figure 4. (Color online). Smoothed profile of the frozen mode in periodic layered 
structures composed of different number TV of unit cells L. The frequency is equal 
to that of the degenerate band edge. The initial rate of growth of the frozen mode 
amplitude is virtually independent of N. The limiting case (d) of the semi-infinite 
structure is also shown in Fig. |3jc). In all cases, the incident wave has the same 
polarization and unity amplitude. The distance z from the surface is expressed 
in units of L. 
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Figure 5. (Color online) Steady-state field distribution in a lossless periodic finite 
slab with Af = 32 at a series of equidistant frequencies lying between the DBE 
frequency oj^ and the frequency u!r of the giant slow wave resonance. Observe 
that the field amplitude increases dramatically as the frequency changes from that 
of the DBE (the lowest solid blue curve) to that of the giant slow wave resonance 
(the upper solid red curve). The lowest solid curve is identical to that in Fig. 

Hb). 
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Figure 6. (Color online). Destructive interference of the propagating and 
evanescent components of the transmitted wave inside semi-infinite photonic 
crystal. The frequency is close but not equal to that of stationary inflection point, 
(a) The squared modulus of the resulting transmitted field - its amplitude at 2 = 
is small enough to satisfy the boundary conditions; (b) the squared modulus of 
the propagating contribution, which is independent of z; (c) the squared modulus 
of the evanescent contribution, which decays with the distance z. The amplitude 
of the incident wave is unity. The distance z from the surface is expressed in units 
of L. 




Figure 7. (Color online). Destructive interference of the two Bloch components 
of the transmitted wave inside semi-infinite photonic crystal. The frequency is 
a; = Wd — lQ~^c/L, which is slightly below the degenerate band edge in Fig. [He), 
(a) The squared modulus of the resulting transmitted field - its amplitude at 
2: = is small enough to satisfy the boundary conditions II17I I: (b) the squared 
modulus of the propagating contribution, which is independent of z; (c) the 
squared modulus of the evanescent contribution, which decays with the distance 
2. The amplitude of the incident wave is unity. 
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Figure 8. (Color online). Destructive interference of the two Bloch components 
of the transmitted wave inside semi-infinite photonic crystal. The frequency is 
oj = ujd + 10~^c/L, which is just above the degenerate band edge, (a) The 
squared modulus of the resulting transmitted field - its amplitude at 2 = is 
small enough to satisfy the boundary conditions (1171 : (b) and (c) the squared 
moduli of the two evanescent contributions; both decay with the distance z. The 
amplitude of the incident wave is unity. 



